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Abstract 

A derivation of stochastic Schrodinger equations is given using quantum filtering theory. We 
study an open system in contact with its environment, the electromagnetic field. Continuous 
observation of the field yields information on the system: it is possible to keep track in real time 
of the best estimate of the system's quantum state given the observations made. This estimate 
satisfies a stochastic Schrodinger equation, which can be derived from the quantum stochastic 
differential equation for the interaction picture evolution of system and field together. Throughout 
the paper we focus on the basic example of resonance fluorescence. 



1 Introduction 

It has long been recognized that continuous time measurements can not be described by the standard 
projection postulate of quantum mechanics. In the late 60's, beginning 70's, Davies developed a theory 
for continuous time measurement culminating in his book [14j . His mathematical work became 
known to the quantum optics community through the paper with Srinivas on photon counting (51 j . 

The late 80's brought renewed interest to the theory of continuous time measurement. For instance the 
waiting time distribution of fluorescence photons of a two-level atom driven by a laser was obtained 
by associating a continuous evolution to the atom in between photon detections and jumps at the 
moments a photon is detected In this way every record of photon detection times determines 

a trajectory in the state space of the atom. Averaging over all possible detection records leads to 
the well-known description of the dissipative evolution of the atom by a master equation. Advantage 
of the trajectory approach is the fact that an initially pure state will remain pure along the whole 
trajectory. This allows for the use of state vectors instead of density matrices, significantly speeding 
up computer simulations CHI; M- 

Infinitesimally, the quantum trajectories are solutions of a stochastic differential equation with the 
measurement process as the noise term. The change in the state is given by the sum of two terms: a 
deterministic one proportional with dt and a stochastic one proportional to the number of detected 
photons dNt in the interval dt. For other schemes such as homodyne detection the corresponding 
stochastic differential equation is obtained as the diffusive limit of photon counting where the jumps 



in the state space decrease in size but become increasingly frequent [HIj this limit the 

stochastic term in the differential equation is replaced by a process with continuous paths. 

The stochastic Schrodinger equations obtained in this way had been postulated before by Gisin |17| . 
|18| . jl5| . in an attempt to generalize the customary unitary evolution in quantum mechanics. The 
stochastic terms are seen as randomness originating from the measurement process. However, in 
this approach the correspondence between the different quantum state diffusion equations and the 
measurements that can be performed is not emphasized. 

Another approach originated from the development of quantum stochastic calculus |201 , 121] , generali- 
zing the classical Ito table to quantum noises represented by creation and annihilation operators (see 
Section 13 . Barchielli saw the relevance of this new calculus for quantum optics 1 . Indeed, in the 
Markovian approximation the interaction between a quantum system and the electromagnetic field is 
governed by a unitary solution of a quantum stochastic differential equation in the sense of |2(J| . 

Belavkin was the first to see the connection between quantum measurement theory and classical 
filtering theory 22 , in which one estimates a signal or system process when observing a function of 
the signal in the presence of noise. This is done by deriving the filtering equation which is a stochastic 
differential equation for the expectation value of the system process conditioned on outcomes of the 
observation process. Belavkin extended the filtering theory to allow for the quantum noises of 
|20| . Stochastic Schrodinger equations turn out to be examples of the quantum filtering or Belavkin 
equation 0, [7]. 

Aim of this paper is to give an elementary presentation of quantum filtering theory. We construct the 
expectation of an observable conditioned on outcomes of a given measurement process. The differential 
form of this conditional expectation is the stochastic Schrodinger equation associated with the given 
measurement. At the heart of the derivation lies the Ito table of quantum stochastic calculus enabling 
a fast computation of the equation. The procedure is summarized in a small recipe in Section [3 

To illustrate the theory we consequently focus on the basic example of resonance fiuorescence of a 
two-level atom for which we consider photon counting and homodyne detection measurement schemes. 
The stochastic Schrodinger equations for these examples are derived in two ways, once via the usual 
approach using quantum trajectories and a diffusive limit, and once using quantum filtering theory. 
In this way we hope to emphasize how conceptually different both methods are. 

This paper is organised as follows. Sections |21 and 13 serve as an introduction to the guiding example 
of this paper: resonance fluorescence of a two-level atom driven by a laser. In Section |2] we put 
the photon counting description of resonance fluorescence by Davies [Hj, 52 1 CHI into the form of a 
stochastic differential equation driven by the counting process. In Section|21we discuss the homodyne 
detection scheme as a diffusive limit of the photon counting measurement, arriving at a stochastic 
differential equation driven by a diffusion process. The equations of Sections |21 and |21 will be rederived 
later in a more general way using quantum filtering theory. 

In Section 0] we introduce the concept of conditional expectation in quantum mechanics by first 
illustrating it in some simple, motivating examples. Section^ldescribes the dissipative evolution of the 
open system within the Markov approximation. The joint evolution of the system and its environment, 
the quantized electromagnetic field, is given by unitaries satisfying a quantum stochastic differential 
equation. Given a measurement of some field observables it is shown how to condition the state of 
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the system on outcomes of the measurement using the construction of Section^ Section |^ is a short 
review of quantum stochastic calculus and its applications to open systems. Sections |51 and El describe 
dilation theory and quantum stochastic calculus in a nutshell. 

Section[71contains the derivation of the quantum filtering equation, the stochastic differential equation 
for the conditional expectation. This equation is the stochastic Schrodinger equation for the given 
measurement. This part ends with a recipe for computing stochastic Schrodinger equations for a large 
class of quantum systems and measurements. The end of the article connects to Sections |21 and El by 
showing how the recipe works in our main example. 



2 The Davies process 



We consider a two-level atom in interaction with the quantized electromagnetic field. The state of 
the atom is described by a 2 x 2-density matrix p, i.e. p > 0, and Trp = 1. Atom and field together 
perform a unitary, thus reversible evolution, but by taking a partial trace over the electromagnetic 
field we are left with an irreversible, dissipative evolution of the atom alone. In the so called Markov 
limit it is given by a norm continuous semigroup {Tt}t>o of completely positive maps. A central 
example discussed in this paper is resonance fluorescence. Here the atom is driven by a laser on the 
forward channel, while in the side channel a photon counting measurement is performed. For the time 
being we will suppress the oscillations of the laser for reasons of simplicity. In this case the Lindblad 
generator of Tt, or Liouvillian L is given by (cf. 9 ): 



d_ 

dt 



t=o 



Tt{p) = L{p) = -i[H,p] + i-[V + V\p]--{V*V,p} + VpV\ where V^r[ ^1, (2.1) 



H := ^cTz is the Hamiltonian of the atom, and is the Rabi frequency. 

The master equation (|2.1() can be unravelled in many ways depending on what photon detection 
measurement is performed. By unravelling the master equation we mean writing L as the sum C + J^, 
where J represents the instantaneous state change taking place when detecting a photon, and £ 
describes the smooth state variation in between these instants. The unravelling for photon counting 
in the side channel is given by (51 

Cip) = ~i[H,p]+i^[V + V*,p] - ^{V*V,p} + (1 - \K,\^)VpV* and Jip) = 

with |ksP the decay rate into the side channel. 

An outcome of the measurement over an arbitrary finite time interval [0, t) is the set of times 
{ti,t2, ■ ■ ■ ,tk} at which photons are detected in the side channel of the field. The number of de- 
tected photons can be arbitrary, thus the space of outcomes is 

oc oo 

n([0,t)):= |Ja„([0,t))= |J{ac[0,i); \a\ ^ n} 

71=0 ri=0 

also called the Guichardet space jl9| . In order to describe the probability distribution of the outcomes 
we need to make ([0,t)) into a measure space. Let us consider the space of n-tuples [0,f)" with its 
Borel (T- algebra and the measure -^Xn where A„ is the Lebesgue measure. Then the map 

jn : [o,i)" 3 (ti,...,t„) {ti,...,t„} e n„{[o,t)) 
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induces the cr-algebra S„ {[0,t)) and the measure /i„ on f7„ ([0,t)). We define now the measure fj, on 
f)([0,t)) such that m({0}) = 1 and /i = /i„ on f7„ ([0,t)). We will abbreviate ^{[0,t)) and S([0,t)) to 
ri* and S*, respectively. 

Davies was the first to show |I31 (see also 101,13) that the unnormalized state of the 2-level atom at 
time t with initial state p, and conditioned on the outcome of the experiment being in a set _E G S* 
is given by: 

M'[E]{p) - / Wt{Lu)ip)dpiu;), 



where ior uj = {ti, . . . ,tk} E with < ti < . . . < tk < t we have 

Wt{Lu)ip) := exp ((t ~tk)C)j.. . J exp ((iz - ti)C)jexp {tiC){p). 

Furthermore, Pp[i?] '■— Tr(A^*[i?](p)) is the probability that the event E occurs if the initial state is 
p. The family of prabability measures {Pp}t>o is consistent, i.e. P*+''[£;] = Pp[i?] for aU € E*, s > 0, 
see hence by Kolmogorov's extension theorem it extends to a single probability measure Pp on the 
cr-algebra of the set 17°°. 

On the measure space (ri°°, Pp) we define the following random variables: 

Nt: 0.°° ^n-. u^\uJC^ [0,t)|, 

counting the number of photons detected in the side channel up to time t. The counting process 
{Nt]t>o has differential dNt := Nt+dt — Nt satisfying dNt{uj) = 1 ii t E uj and dNt{uj) — otherwise. 
Therefore we have the following Ito rules: dNtdNt — dNt and dNfdt = 0, (cf. .2.). 
To emphasise the fact that the evolution of the 2-level atom is stochastic, we will regard the normalized 
density matrix as a random variable {p*t\t>o with values in the 2 x 2-density matrices defined as follows: 

p\ : 17- ^ Af2 : ^ - := /\ I ' ^M^^ (2.2) 

Tr(m(wn[0,t))(p)j 

The processes Nt and p\ are related through the stochastic differential equation dp\ ~ atdt + PtdNt- 
Following [2] we will now determine the processes at and (3t by differentiating (|2.2|l . If < e then 
dNt{uj) = 1, i.e. the differential dt is negligible compared to dNt — 1, therefore: 

Tr(J^(pL)) 

On the other hand, \i t ^ lo then dNt{ijo) — 0, i.e. dNt is negligible compared to dt. Therefore it is 
only the dt term that contributes: 



, . d exp{{s~t)C){p'J 



ds 



'Tr(^exp{{s-t)C){pl) 



^fl) - ^r^Tr(£(pL)) = C{pl) + Tr(^(pL))pL, 



(2.4) 



where we used that Tr(£(p^)) ^Tr(j'(/9^)) , as a consequence of the fact that Tr(L(cr)) = for 
all density matrices a. Substituting 1)2.3(1 and (|2.4|l into dp^ = atdt + PtdNt we get the following 
stochastic Schrodinger equation for the state evolution of the 2-level atom if we are counting photons 
in the side channel (cf. |ll)p: 



dpi = L{pi)dt + (^r^^ - pi) {dNt - Tr{Jipi))dt 



(2.5) 
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The differential dMt :— dNt — Tr(j'(p* and the initial condition Mq — define an important 
process Mt called the innovating martingale, discussed in more detail in Sectional 



3 Homodyne detection 

We change the experimental setup described in the previous section by introducing a local oscillator, 
i.e. a one mode oscillator in a coherent state given by the normalised vector in P{N) 

^(a,) := exp (^^) (1, a*, ^, . . .), (3.1) 

for a certain at € C. We take at = — , where Wt is a complex number with modulus \wt\ — 1, and 
e > 0. The number e is inversely proportional to the intensity of the oscillator. Later on we will let 
the intensity go to infinity, i.e. £ ^ 0. The phase (pt of the oscillator is represented by Wt = exp(z0t), 
with (f)t = (j)o + ujiot, where ujio is the frequency of the oscillator. 

The local oscillator is coupled to a channel in the electromagnetic field, the local oscillator beam. 
The field is initially in the vacuum state. The local oscillator and the field are coupled in such a way 
that every time a photon is detected in the beam, a jump on the local oscillator occurs, given by the 
operation 

Jioip) = AiopAl, (3.2) 

where Aio is the annihilation operator corresponding to the mode of the local oscillator. The coherent 
state ipi^tt) is an eigenstate of the jump operator Aio at eigenvalue at- 

Now we are ready to discuss the homodyne detection scheme. Instead of directly counting photons 
in the side channel we first mix the side channel with the local oscillator beam with the help of a 
fifty-fifty beam splitter. In one of the emerging beams a photon counting measurement is performed. 
A detected photon can come from the atom through the side channel or from the local oscillator via 
the local oscillator beam. Therefore the jump operator on states a of the atom and the oscillator 
together, is the sum of the respective jump operators: 

Ja®io(o-) = [KsV ® I + I ® Aio)a{KsV* ^ I + I ^ A^J. 

An initial product state p (g) \ip{at)){il'{c(t)\ of the 2-level atom and the local oscillator will remain a 
product after the jump since ijjiat) is an eigenvector of the annihilation operator. Tracing out the 
local oscillator yields the following jump operation for the atom in the homodyne setup: 

Jaip) = Trio(^Ja^io{p® |V'(at))(^(at)|)) = {k,V + ^) p{ksV* + ^). 

In the same way as in Sectional we can derive the following stochastic Schrodinger equation for the 
state evolution of the two-level atom when counting photons after mixing the side channel and the 
local oscillator beam [2] |1U| : 

dpi = L{pi)dt+^[,^0^ pi)e[dNt - Tr{Mpi))dt), (3.3) 

where the extra e's are introduced for future convenience. We will again use the abbreviation: dMf = 
dNt — Tr(j7a(/5i))di for the innovating martingale (see Sectiond). In the homodyne detection scheme 
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the intensity of the local oscillator beam is taken extremely large, i.e. we are interested in the limit 
e ^ 12], 0, ESI- Then the number of detected photons becomes very large and it makes sense 
to scale and center Nt, obtaining in this way the process with differential dW^ := edNt — dt/e and 
Wq — 0. We find the following Ito rules for dW^: 

dW^dW^ = {edNt ~ ^dt) {edNt - jdt) = e^dNt = edW^ + dt, 
dW^dt ^ 0. 

In the limit £ — > this becomes dWtdWt = dt and dWtdt ~ 0, i.e. the process Wt '■— linij^o Wt is a 
diffusion. It is actually this scaled and centered process that is being observed and not the individual 
photon counts Nt, see We pass now to the evaluation of the limit of (|3.3() : 

lim -( f;^^'^ - pi) = wtK.piV* +Wtn,Vpi - Tr{wtT^,plV* +WtK,Vpl)pl. 

This leads to the following stochastic Schrodinger equation for the homodyne detection scheme 

uni, m 

dpi = L{pl)dt + {wtlisplV* + wt^isVpi - Tr{wtKspiV* + WtKsVpDpl ) dAf/'^ (3.4) 
for all states p G M2, where 

dMf ^ dWt - TT{wtK,plV* + WtKsVpl)dt. (3.5) 



Let as{t) and af,{t) denote the annihilation operators for the side channel and the local oscillator 
beam, respectively. They satisfy the canonical commutation relations 

[ai{t),a*{r)] = 6i^jS{t~r), i,j G {s,b}. 

Smearing with a quadratically integrable function / gives 

A{f) = J .f{t)adt)dt, ie{s,b}. 

By definition, the stochastic process {Nt}t>o counting the number of detected photons has the same 
law as the the number operator A{t) up to time t for the beam on which the measurement is performed. 
Formally we can write 

A{t)= I {a*{r)(E)I + I(g)al{r)){as{r)(g)I + I(g)ab{r))dr. 
Jo 

The oscillator beam is at time t in the coherent state "0 ("7)1 where ft € L^(R) is the function 

r i"> ifrX[o.t] (^)- Since the state of the local oscillator beam is an eigenvector of the annihilation 
operator afc(r) 

, s , / /* \ Wr ( ft 



^e ) e 
we find 

em --= eA^t) ®I + e f (— ^(r) + ^a,(r)) ® / + - - 

e Jq £ £ £^ £ 

= eA,{t) ®I+ {a: (ft) + A.ift)) ® /. 

The operator X^{t) :— A*(/t) + As{ft) is called a field quadrature. We conclude that in the limit 
£ — > the homodyne detection is a setup for continuous time measurement of the field quadratures 
X^{t) of the side channel, (cf. *^). 
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4 Conditional expectations 



In the remainder of this article we will derive the equations H2.5|l and (|3.4II in a different way. We 
will develop a general way to derive Belavkin equations (or stochastic Schrodinger equations). The 
counting experiment and the homodyne detection experiment, described in the previous sections, serve 
as examples in this general framework. The method we describe here closely follows Belavkin's original 
paper on quantum filtering theory |^. Our approach differs in its construction of the conditional 
expectation, which is the topic of this section. 

Let us remind the concept of conditional expectation from probability theory. Let S,P) be a 
probability space describing the "world" and E' C E a a-algebra of events to which "we have access" . 
A random variable / on (fi, E,P) with E(|/|) < oo can be projected to its conditional expectation 
£{f) which is measurable with respect to E' and satisfies 

/ /dP= / £{f)dP 

J E J E 

for all events S in E'. Our information about the state of that part of the world to which we have 
access, can be summarized in a probability distribution Q on E'. Then the predicted expectation of / 
given this information is £'(/)dQ. We will extend this now to quantum systems and measurements. 

The guiding example is that of an n level atom described by the algebra B := M„ undergoing a 
transformation given by a completely positive unit preserving map T : B B with the following Kraus 
decomposition T{X) — '^i^fiV*XVi. The elements of $7 can be seen as the possible measurement 
outcomes. For any initial state p oi B and measurement result i € n, the state after the measurement 
is given by 

= V,pV*/TT{V,pVn, 

and the probability distribution of the outcomes is p = X^iGnP'^* where Si is the atomic measure at 
z, and Pi — TT{VipVi*), which without loss of generality can be assumed to be strictly positive. We 
represent the measurement by an instrument, that is the completely positive map with the following 
action on states 

M : M* ^ M* ®£\n): p^Y,p,(g> pA. (4.1) 

Let X G S be an observable of the system. Its expectation after the measurement, given that the 
result i £ D, has been obtained is Tv{piX). The function 

£{X) : n^C: Tr(p,X) 

is the conditional expectation of X onto £°°{fl). li q ~ "^QiSi is a probability distribution on 17 then 
'^qi£{X){i) represents the expectation of X on a statistical ensemble for which the distribution of 
the measurement outcomes is q. We extend the conditional expectation to the linear map 

£ : B(g)e°°{n) £°°{n) c B(g)£°^{n) 

such that for any element A : i i-^ Ai in B ^ £°°{n) ^ £°°{n ^ B) we have 

£{A): Tr(p,A,). 

This map has the following obvious properties: it is idempotent and has norm one. Moreover, it is 
the unique linear map with these properties preserving the state Ai{p) on B®£°°{Vl). For this reason 
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we will call £, the conditional expectation with respect to A4{p). Its dual can be seen as an extension 
of probability distribitions q £ £^{n) to states on B (E) £°°{fl) 



Thus while the measurement H4.1|l provides a state Ai{p) on B iS) the conditional expectation 

with respect to Ai{p) extends probability distributions q G £^{^l) of outcomes, to states on B ® 
and in particular on B which represents the state after the measurement given the outcomes 
distribution q. 

With this example in mind we pass to a more general setup which will be needed in deriving the 
stochastic Schrodinger equations. Let ^ be a unital *-algebra of bounded operators on a Hilbert 
space 11 whose selfadjoint elements represent the observables of a quantum system. It is natural 
from the physical point of view to assume that A is strongly closed, i.e. if {j4„}„>o is a sequence of 
operators in A such that HAn?/)!! 11^-011 for any vector -0 in H and a fixed bounded operator A, then 
A £ A. From the mathematical point of view this leads to the rich theory of von Neumann algebras 
inspired initially by quantum mechanics, but can as well be seen as the generalization of probability 
theory to the non-commutative world of quantum mechanics. Indeed, the building blocks of quantum 
systems are matrix algebras, while probability spaces can be encoded into their commutative algebra 
of bounded random variables L°°{il, S,P) which appeared already in the example above. A state is 
described by a density matrix in the first case or a probability distribution in the second, in general 
it is a positive normalized linear functional : A C which is continuous with respect to the 
weak*-topology, the natural topolgy on a von Neumann algebra seen as the dual of a Banach space 



Definition 4.1: Let B he a, von Neumann subalgebra of a von Neumann algebra A of operators 
on a (separable) Hilbert space H. A conditional expectation of A onto S is a linear surjective map 
£ : A^ B, such that: 

1. £^ — £ {£ is idempotent), 

2. Waga ■ < \\M {£ is normcontractive). 



In it has been shown that the conditions 1 and 2 are equivalent to £ being an identity preserving, 
completely positive map, and satisfying the module property 



generalizing a similar property of conditional expectations in classical probability theory (cf. 

In analogy to the classical case we are particularly interested in the conditional expectation which 
leaves a given state p on A invariant, i.e. po £ = p. However such a map does not always exist, but if 
it exists then it is unique [H2I and will be denoted £p. Using £p we can extend states u on S to states 
a o £p oi A which should be interpreted as the updated state of A after receiving the information (for 
instance through a measurement) that the subsystem B is in the state a (cf. j25i|). 
In the remainder of this section we will construct the conditional expectation £p from a von Neumann 
algebra A onto its center C := {C E A; AC = CA for all A £ A} leaving a given state p on A 




£{BiAB2) = Bi£{A)B2, 



for all Bi,B2 e B, and A e A, 



(4.2) 
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invariant. The center C is a commutative von Neumann algebra and is therefore isomorphic to some 
L°°{^l,I],P). In our guiding example the center of B^£°°{n) is £°°{^l). Later on (see section^ this 
role will be played by the commutative algebra of the observed process with il the space of all paths 
of measurement records. 



Theorem 4.2: There exists a unique conditional expectation £p : C which leaves the state p on 
A invariant. 

Proof. The proof is based on the central decomposition of A In our guiding example, BiS)i°°{^) 
is isomorphic to (Bi^nBi where the Bi^s are copies of B. In general we can identify the center C with 
some L°°{^1, S,P) where P corresponds to the restriction of p to C. We will ignore for simplicity all 
issues related with measurability in the following constructions. The Hilbert space H has a direct 
integral representation H = ]HIi^P(dw) in the sense that there exists a family of Hilbert spaces 
{Ht^jijgo and for any e H there exists a map oj h-s- ^/^^^ g such that 

i^A) = / (V^,(/.^)P(du;). 
Jn 

The von Neumann algebra A has a central decomposition A — ALuPidoj) in the sense that there 
exists a family {^c^jc^gn of von Neumann algebras with trivial center, or factors, and for any A € A 
there is a map oj i-^ A^j G Au; such that {Aip)uj — A^^ipi^ for all S H and P-almost all u! G ft. The 
state p on A has a decomposition in states p^j on A^j such that for any A e ^ its expectation is 
obtained by integrating with respect to P the expectations of its components A^^ : 

p{A) = ( p^{A^)¥{dw). (4.3) 
Jn 

The map £p : A ^ C defined by 

£p{A) : uj Pu^iAu;) 

for all A €z A is the desired conditional expectation. One can easily verify that this map is linear, 
identity preserving, completely positive (as a positive map onto a commutative von Neumann algebra), 
and has the module property. Thus, £p is a conditional expectation and leaves the state p invariant 
bv l4.3l Uniqueness follows from [321 • □ 

It is helpful to think of the state p and an arbitrary operator A as maps p, : u t—^ puj, and respectively 
At : UJ i—f A^. The conditional expectation £p{A) is the function (A,) : lu t-^ Pu{Ai^). 



5 The dilation 

Let B be the observable algebra of a given quantum system on the Hilbert space H. In the case of 
resonance fluorescence B will be all 2 x 2 matrices M2, the algebra of observables for the 2-level atom. 
The irreversible evolution of the system in the Heisenberg picture is given by the norm continuous 
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semigroup {Tt}t>o of completely positive maps Tt : B ^ B. By Lindblad's theorem [221 we have 
Tt — exp{tL) where the generator L : B ~> B has the following action 

k 

L{X) = i[H,X]+Y.^;XV, - -{V*V,,X}, (5.1) 
where H and the V,'s are fixed elements of S, H being selfadjoint. 

We can see the irreversible evolution as stemming from a reversible evolution of the system B coupled 
to an environment, which will be the electromagnetic field. We model a channel in the field by the 
bosonic or symmetric Fock space over the Hilbert space L^(R) of square integrable wave functions on 
the real line, i.e. 

oo 

The algebra generated by the field observables on T contains all bounded operators and we denote it 
by W. For the dilation we will need k independent copies of this algebra W®'^. 

The free evolution of the field is given by the unitary group St, the second quantization of the left 
shift s{t) on L^(]R) , i.e. s{t) '■ f ^ f{- + 1). In the Heisenberg picture the evolution on W is 

SlWSt -.^ AA[St]{W). 

The atom and field together form a closed quantum system, thus their joint evolution is given by a 
one-parameter group {Tt^tes. of ^-automorphisms on B ® yV"^*=: 

X ^ U;XUt := Ad[Ut]{X). 

The group Ut is a perturbation of the free evolution without interaction. We describe this perturbation 
by the family of unitaries Ut :— S^t^t for alH S R satisfying the cocycle identity 

Ut^^ ^ S'!!^UtSf''Us, for all t,se«.. 

The direct connection between the reduced evolution of the atom given by 15. 1|) and the cocycle Ut 
is one of the important results of quantum stochastic calculus '20' which makes the object of Section 
|S1 For the moment we only mention that in the Markov limit, Ut is the solution of the stochastic 
differential equation PU]. f ^ . 

dUt = {V,dA*{t) - V*dA,it) - (iH + ^V*Vj)dt}Ut, Uo = 1, (5.2) 

where the repeated index j is meant to be summed over. The quantum Markov dilation can be 
summarized by the following diagram (see |24j ) : 

B B 

Id®i«"=| |ld«i0®'= (5.3) 

S^yy®fe > B(»yV^'' 

i.e. for all X £ B : Tt{X) ^ (id (g) cj)'^'') {ft(X (g) l'^'^)), where (/) is the vacuum state on W, and 1 is 
the identity operator in W. Any dilation of the semigroup Tt with Bose fields is unitarily equivalent 
with the above one under certain minimality requirements. The diagram can also be read in the 



10 



Schrodinger picture if we reverse the arrows: start with a state p of the system B in the upper right 
hand corner, then this state undergoes the following sequence of maps: 

p p ^ (p o ft = ft,{p ® </)®^-) ^ Tr^«. {ft, [p ® 0®'=)) . 

This means that at < = 0, the atom in state p is coupled to the k channels in the vacuum state, and 
after t seconds of unitary evolution we take the partial trace taken over the k channels. 

We would now like to introduce the measurement process. It turns out that this can be best described 
in the interaction picture, where we let the shift part of Ut — Sf^Ut act on the observables while the 
cocycle part acts on the states: 

p\X):^p®4>®^{U:XUt) (5.4) 

for &\\ X ^ B ® yV"**^ . It is well known that for the Bose field for arbitrary time t we can split the 
noise algebra as a tensor product 

w = Wq) ® yV[o,t) ® W[t 

with each term being the algebra generated by those fields over test functions with support in the 
corresponding subspace of L^(R): 

L2(R) = L^{{-oo,Q))®L^{[Q,t))®L^{[t,oo)). 

Such a continous tensor product structure is called a filtration and it is essential in the development 
of quantum stochastic calculus reviewed in Section El The observables which we measure in an 
arbitrary time interval [0, t) form a commuting family of selfadjoint operators {Ys\Q<s<t whose spectral 
projections belong to the middle part of the tensor product W[o,t)- In the Davies process Ys — A(s), 
i.e. the number operator up to time s, while in the homodyne case Yg = X^{s). Notice that the 
part Wg) will not play any significant role as it corresponds to "what happened before we started our 
experiment" . 

Let Ct be the commutative von Neumann generated by the observed process up to time {l's}o<s<t 
(t > 0), seen as a subalgebra oiB® W"**^. By a theorem on von Neumann algebras, Ct is equal to the 
double commutant of the observed process up to time t: Ct — {Ys\ < s < t}", with the commutant 
S' of a subset 5 of 6 « W^^ being defined by 5' := {X € 6 (g) W®'=; XS ^ SX yS € S}. The 
algebras {Ct}t>o form a growing family, that is Cg C Ct for all s < t. Thus we can define the inductive 
limit Coo := limt^oo Ct, which is the smallest von Neumann algebra containing all Ct- On the other 
hand for each i > we have a state on Ct given by the restriction of the state p* of the whole system 
defined by (|5.4(l . We will show now that that the states p* for different times "agree with each other" . 

Theorem 5.1: On the commutative algebra Coo there exists a unique state p°° which coincides with 
p* when restricted to Ct C Coo, for all t>Q. In particular there exists a measure space (fi, S, Pp) such 
that (Coo,p°°) is isomorphic with L°°{n,'E,Fp) and a growing family {T,t}t>o of a-subalgebras of S 
such that {Ct,p^) = L°°(f7,St,Pp). 



Proof. In the following we will drop the extensive notation of tensoring identity operators when 
representing operators in W[s,t) for all s,i e M. Let X e Cg, in particular X € Wj^^-,. By (|5.2|l . 

Ut € B 1^ W^\y because the coefficients of the stochastic differential equation lie in B ^[o *)■ This 

implies that S'^^UtSf' e B®W^^^^^^y Using the tensor product structure of W®^ we see that W^'^^^ 



11 



and B (g) j,^^-, commute, and in particular X commutes with S^^UtSf^. Then 

p'+'{x) = p^iu^^^xut+s)^ p''iu:is'^'^UtSf''rxs^^UtSf''Us) 

= p%U:XUs)^p'{X). (5.5) 

This imphes that the hmit state p°° on Coo with the desired properties exists, in analogy to the 
Kolmogorov extension theorem for probability measures. As seen in the previous section, (Coo,p°°) 
is isomorphic to L°°(ri, E,Pp) for some probabihty space (fi, E,Pp). The subalgebras (Ct,p*) are 
isomorphic to L°°(17, Ej, Pp) for some growing family {Et}t>o of cr-subalgebras of E. □ 

Remark. From spectral theory it follows that the measure space (ri*,Ej) coincides with the joint 
spectrum of {Ys}s<t, i-e. il* is the set of all paths of the process up to time t. For the example of the 
counting process this means that Jl* is the Guichardet space of the interval [0,t), which is the set of 
all sets of instants representing a "click" of the photon counter, i.e. it is the set of all paths of the 
counting process. 

We define now At ■= C[ for all t > , i.e. At is the commutant of Ct, then Ct is the center of the 
von Neumann algebra At ■ Notice that the observable algebra of the atom B is contained in At ■ By 
Theorem 14.21 we can construct a family of conditional expectations ; At — > Ct}t>t)- For each 
t, f*t depends on the state of the "world" at that moment p*, keeping this in mind we will simply 
denote it by £*. An important property of is that p°° o £* = p^ o £* = p*^ since the range of is 
Ct and leaves invariant. 

For an element X e At, £^{X) is an element in Ct, i.e. a function on fit. Its value in a point lo G fit, 
i.e. an outcomes record up to time t, is the expectation value of X given the observed path lj after t 
time units. We will use the notation £*'{X) :— (X,) defined in the end of Section^jto emphasise the 
fact that this is a function on fif . When restricted to B ®Ct the conditional expectation is precisely 
of the type discussed in our guiding example in Section^ 

There exists no conditional expectation from B ®W onto Ct since performing the measurement has 
demolished the information about observables that do not commute with the observed process jS]. 
We call At the algebra of observables that are not demolished ^ by observing the process {Kj}o<s<t- 
This means that performing the experiment and ignoring the outcomes gives the same time evolution 
on At as when no measurement was done. 

From classical probability it follows that for all t > there exists a unique conditional expectation 
lEp : Coo Ct that leaves the state p°° invariant, i.e. p°° o Ep = p°° . These conditional expectations 
have the tower property, i.e. Ep o Ep = Ep for all t > s > 0, which is often very useful in calculations. 
Ep is the expectation with respect to Pp, and will simply be denoted Ep. Note that the tower property 
for s = is exactly the invariance of the state p°°{— Ep). 



6 Quantum stochastic calculus 

In this section we briefly discuss the quantum stochastic calculus developed by Hudson and Parthasarathy 
|2U| . For a detailed treatment of the subject we refer to ^9J and Let !F{Ti,) denote the symmetric 
(or bosonic) Fock space over the one particle space Ti := 'C^ ® i^(]R+) = i^({l, 2, . . . fc} x M+). The 
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space C*^ describes the k channels we identified in the electromagnetic field. As in the previous section 
we denote the algebra of bounded operators on the one channel Fock space J^(]R+) by W, and on the 
k channels T{n) by 

For every / e 7i we define the exponential vector e(/) G J'i'H) in the following way: 

e(/) :=l®0^/«", 

n—l ^ 

which differs from the coherent vector by a normalization factor. The inner products of two exponential 
vectors e(/) and e{g) is {e{f),e{g)) = exp((/, g)). Note that the span of all exponential vectors, 
denoted T>, forms a dense subspace of J^{7i). Let fj he the j'th component of / G H for j = 1, 2, . . . , fc. 
The annihilation operator Aj(t), creation operator A*(t) and number operator Aij{t) are defined on 
the domain V by 

MMf) = {X[o,t], fMf) = f fiis)ds e{f) 

Jo 

{e{g), A*{t)e{f)) = {g„ X[o.t]>(e(5), - f g,is)ds exp((/,g)) 

Jo 

(e(g), Ay(i)e(/)) = {g„ X[o.t]fj){e{g), e(/)) = / g,{s)fj{s)ds exp((/,g)). 

Jo 

The operator Aii{t) is the usual counting operator for the i'th channel. Let us write L^(R+) as direct 
sum i2([0,t]) e L2([i,oo]), then J^{L'^{R+)) is unitarily equivalent with T{L'^{[0,t]) (8> J^{L'^[t, oo)) 
through the identification e(/) = e(/t]) (g) e(/[(), with /«] = /x[o,t] and f[t = /x[t,oo)- We will also 
use the notation f[s.t] for fx[s.t] and omit the tensor product signs between exponential vectors. The 
same procedure can be carried out for all the k channels. 

Let Mf be one of the processes Aj{t),A*{t) or A,j{t). The following factorisability property [201, EHl 
makes the definition of stochastic integration against Mj possible 

{Mt - Ms)e{f) = e(/,]){(Afi - M,)e(/[,,,])}e(/[,), 

with {Alt — AIs)e{f[s.t]) & ^{C' (8" L'^{[s, t])) . We firstly define the stochastic integral for the so called 
simple operator processes with values in the atom and noise algebra B (g> >V®'= where B := M„ 

Definition 6.1: Let {Ls}o<s<t be an adapted (i.e. Lg G B^Ws] for all < s < simple process with 
respect to the partition {sq = 0, Si, . . . , Sp — t} in the sense that Lg ~ Lg. whenever Sj < s < sj^i. 
Then the stochastic integral of L with respect to AI on C" ® I? is given by [20] , [201 ■ 

t p-i 
LsdAI, fe{u) ^ {L,Je{us^])){{AI,^^, - M,Je(w[,^,,^^,]))e(M[,^^ J. 



By the usual approximation by simple processes we can extend the definition of the stochastic integral 
to a large class of stochastically integrable processes [201; [20]' We simplify our notation by writing 
dXt = LtdAIt for Xt = Xq + LsdAIg. Note that the definition of the stochastic integral implies that 
the increments dMs lie in the future, i.e. dAIg G 'W[s- Another consequence of the definition of the 
stochastic integral is that its expectation with respect to the vacuum state (j) is always due to the 
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fact that the increments dAj,dA*,dAij have zero expectation values in the vacumir. This wiU often 
simphfy calculations of expectations, our strategy being that of trying to bring these increments to 
act on the vacuum state thus eliminating a large number of differentials. 

The following theorem of Hudson and Parthasarathy extends the Ito rule of classical probability theory. 



Theorem 6.2: (Quantum ltd rule '2t)"(, ^29f ) Let Mi and M2 he one of the processes Aj,A* or 
Aij. Then M1M2 is an adapted process satisfying the relation: 

dMiNh = MidAh + MadAfi + dMidlVh, 

where dAIidM2 is given by the quantum ltd table: 



dMi\dM2 


dA* 


dAij 


dAi 


dAl 











dAki 


SlidAl 


SiidAkj 





dAk 


5kidt 


6 kid A j 






Notation. The quantum Ito rule will be used for calculating differentials of products of Ito integrals. 
Let {•^i}i=i,- -,p -^^o integrals, then 

d{ZiZ2...Zj,)= M 

i/C{l,....p} 

where the sum runs over all non-empty subsets of {1, . . . ,p} and for any v = {ii, . . . ik}, the term 
[u] is the contribution to diZxZ^ ■ ■ ■ Zp) coming from differentiating only the terms with indices in 
the set {ii, . . .ik} and preserving the order of the factors in the product. For example the differ- 
ential d{ZiZ2Z3) contains terms of the type [2] = Zi{dZ2)Z3, [13] = {dZi)Z2{dZ3), and [123] = 

idZi)idZ2)idZ3). 

Let Vj for j — 1,2, . . . , k, and H be operators in B with H is selfadjoint. Let 5' be a unitary operator 
on C" (8) /^({l, 2, . . . , fc}) with Sij — (i, Sj) £ B the "matrix elements" in the basis {|i >: i = 1, . . . , fc} 
of C'^. Then there exists a unique unitary solution for the following quantum stochastic differential 
equation [201, EHI 

dUt = {V,dA*{t) + (5y - %)dA,,(t) - V*S,,dA,{t) - {iH + ^V*Vj)dt}Ut, (6.1) 

with initial condition Uq = 1, where again repeated indices have been summed. Equation l|5.2|l . 
providing the cocycle of unitaries perturbing the free evolution of the electromagnetic field is an 
example of such an equation. The terms dAij in equation Hfi.l|l describe direct scattering between the 
channels in the electromagnetic field We have omitted this effect for the sake of simplicity, i.e. we 
always take Sij — Sij . 

We can now check the claim made in Section |5l that the dilation diagram 15.31 commutes. It is easy 
to see that following the lower part of the diagram defines a semigroup on B. We have to show it is 
generated by L. For all X B we have 

dId(8 0'=(ft(X® 1®*^)) =ld(g)(p''{d U;x®l^^Ut). 
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Using the Ito rules we obtain 

d U;X (g) l^'^Ut = {dU*)X (g) I'^'^Ut + U*X (g) l^'^dUt + {dU*)X (g) I'^'^dUt. 

With the aid of the Ito table we can evaluate these terms. We are only interested in the dt-terms since 
the expectation with respect to the vacuum kills the other terms. Then we obtain: d Id (t>^{UiX (g) 
I'^^Ut) = Id ® cj)^[U;L{X) (g> l®''Ut)dt, proving the claim. 

Now we return to the example of resonance fluorescence. Suppose the laser is off, then we have spon- 
taneous decay of the 2-levcl atom into the field which is in the vacuum state. For future convenience 
we already distinguish a forward and a side channel in the field, the Liouvillian is then given by 



L{X)=i[H,X] + v:xv„-\{v:v,,x}, 



where 

1 



with and \ks\'^ the decay rates into the forward and side channel respectively. 

The dilation of the quantum dynamical system (M2, {Tj = exp(tL)}t>o), is now given by the closed 

system (M2 ® W/ (g Wg) {2t}tGR) with unitary cocycle given by 

dUf = {VfdA){t) - V/dAf{t) + VsdA*{t) - V;dAs{t) - {iH + ^V*V)dt}Ut"^, U^'^ = 1, 

where the superscript sd reminds us of the fact that the laser is off, i.e. we are considering spontaneous 
decay. We can summarize this in the following dilation diagram 

^ T,=cxp{tL)^ ^ 



Idi8)i«)i 

Bg>Wfg)Ws Bg)yVf(g,yVs 

where JJ^'^ is given by St <g) StU^'^ for t > 0. 

We change this setting by introducing a laser on the forward channel, i.e. the forward channel is now 
in a coherent state (seel^Jl '■= {"^{h), ■ i^{h)) for some h e L^(R+). This leads to the following 
dilation diagram 



B B 



(6.2) 



B(g>Wf(g)Ws — B(g>Wf(g)Ws 



i.e. the evolution on B has changed and it is in general not a semigroup. Denote by W{h) the unitary 
Weyl or displacement operator defined on V by: W{h)'ip{f) = cxp(— 2ilm(/i, f))'>p{f + h). Note that 
W{h)(j} = W{h)ilj{0) — ip{h), so that we can write 

Tl\X) = Id (g 7,, (g) (t){Uf*X (g 1 (g lU^^) = Id (g g) <j}{Ws{hyU^'^*X g) 1 (g lU^'^Wfih)) = 
Id (g (/) (g (l){Wf{hfO*Uf*X g) 1 (g lUfWf{ht{)), 
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where ht] := hx[o,t] and Wf{h) -.^ 1(g) W(h) g) 1. Defining Ut := U^'^Wf{ht]), together with the 
stochastic differential equation for Wf{hf]) |^ 

dWfiht]) = {h{t)dA}it) - h{t)dAf{t) - ^\h{t)\''dt}W}{hf^, Wfiho) = 1, 

and the Ito rules leads to the following quantum stochastic differential equation for Ut- 
dUt = {{Vf + h{t))dA*{t) ~ {Vf* +h{t))dAf{t) + V,dAl{t) ~ V*dA,{t) - 

-{iu ^^\{\h{t)\^ + v*v + 2h(t)v*))dt )\Jt, c/o = i. 

Define V/ := + V; := K and if := i? + i\(h{t)Vf - h{t)V/) then this reads 

dUt^ J2 {VadAlit)-V;dA,-^iiH + v;K)dt}Ut, Uo^l. (6.3) 

The time dependent generator of the dissipative evolution in the presence of the laser on the forward 
channel is 

L{X) = t[H,X] + J2 t*XV, ~ l{V:Va,X}. (6.4) 
Therefore the diagram for resonance fluorescence (|6.2|) is equivalent to 

B B 

Id®!®! 



Bg)Wf(E)Ws B(E)Wfg)Ws 

where Ut is given by StgjStUt for t > 0. For h{t) = —ifl/nf, we find the master equation for resonance 
fluorescence (I2.1|l . From now on we will no longer suppress the oscillations of the laser, i.e. we take 
h(t) = —ieKp{iLut)Q/Kf. Then we find 

L{X) = i[H,X] - i^[e-"^V + e'"V*,X] - ^{V*V,X} + V*XV, 
note that the laser is resonant when lo — ujq. 



7 Belavkin's stochastic Schrodinger equations 

Now we are ready to derive a stochastic differential equation for the process £*{X). In the next section 
we will see that this equation leads to the stochastic Schrodinger equations (|2.5|) and (|3.4f) . that we 
already encountered in Sections [3 and O 

Definition 7.1: Let X be an element of B Mn- Define the process {M^}t>o in the algebra 
Coc = L°°{il, S],Pp), generated by the observed process {Yt}t>o (see Section|S} by 

M,"" ■.^£\X)-£°{X)~ f £^{L{X))dr, 

Jo 
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where L : B ^ B is the Liouvilhan. In the following we suppress the superscript X in to simplify 
our notation. 



Note that from the above definition it is clear that Mt is an element of Ct for all t > 0. The following 
theorem first appeared (in a more general form and with a different proof) in [S] and is at the heart of 
quantum filtering theory. We prove it using the properties of conditional expectations. For simplicity 
we have restricted to observing a process in the field W®'^. The theory can be extended to processes 
that are in B iS) yV®*^, transforming it into a more interesting filtering theory. For the stochastic 
Schrodinger equations arising in quantum optics our approach is general enough. 

Theorem 7.2: The process {Mt}t>o of deG.nition \7.1\ is a martingale with respect to the Gltration 
{'E,t}t>o of O and the measure Pp, i.e. for all t > s > we have: Ep{Mt) = Mg. 

Proof. From the module property of the conditional expectation it follows that Ep{Mt) — Mg for 
t > s > is equivalent to Wp{Mt — Ms) = for t > s > 0. This means we have to prove for all 
t> s>0 and E eY,^: 

I E^(Mt-M,)HPp(dc^) =0, 
Je 

which, by the tower property, is equivalent to 

/ {Mt~Ms){Lu)Pp{dLo)^0, (7.1) 

JE 

i.e. Ep(x_B(Mf — Ms)) — 0. Now using Definition |0] and again the module property of the conditional 
expectation we find, writing E also for the projection corresponding to xe 

Ep{xE{Mt - Ms)) = *(X (g> E) - £^X (g> E) - £-'{L{X) (g, E)dr^ 

= p\X (g) E) - p'iX (g) E) - p''{L{X)(gE)dr. 

This means we have to prove: dp\X®E)- p\L{X)g)E)dt = 0, for all t > s. Note that p^iXgjE) = 
p°iU^X(gEUt) = p®(l)'^^{U;X®EUt). Therefore dp\X ® E) = p ® (j)®^ {d{U^ X ® EUt)) . We wih 
use the notation below Theorem 16 . 21 with Zi = and Z2 ^ X ® EUt. Using the quantum Ito table 
and the fact that only the dt terms survive after taking a vacuum expectation, we find: 

dp"(t/;X ® EUt) = p"([l]) + P°([2]) + P°([12]), where 

([1]) + ([2]) = (t/; {i[H, X]®E- \{V;V, , X} g> E)Ut) dt 
p°([12]) = p"{u;iV*XV,) g) EUt)dt. 

This means dp\X g)E)= p\L{X) g) E)dt, for all t > s, proving the theorem. □ 

Note that in the proof of the above theorem we have used that the projection E £ Cs commutes 
with the increments dAj(s), dA*(s), ds and with the processes in front of the increments in equation 
(|5.2I) . i.e. V,, V* , V*Vj and H. If the theory is extended to a more general filtering theory 6 , then 



17 



these requirements become real restrictions on the process {Yt}t>Q- If they are satisfied the observed 
process {Yt}t>Q is said to be self non demolition 1^. 

Definition 17.11 impUes the following stochastic differential equation for the process £'^{X) 

d£*{X) ^£\L{X))dt + dMt, (7.2) 

called the Belavkin equation. The only thing that remains to be done is linking the increment dMt to 
the increment of the observed process dYt. 

Let us assume that the observed process {yf}t>o satisfies a quantum stochastic differential equation 
dYt aj{t)dA*{t) + P,j{t)dA,j{t) + a*{t)dAj{t) + 5{t)dt, 

for some adapted stochastically integrable processes aj , , and 5, such that aj {t) , (3ij (t) ,5{t) € W®*^ 
for all t > 0, and p*j = f3ji, 6 = 6* since Yt is selfadjoint. Furthermore, since the observed process 
{Yt}t>o is commutative, we have [dY(,ys] = for all s <t, which leads to 

[a,{t),Ys]dA*it) + [p,,{t),Ys]dA,,{t) + [a*{t),Ys]dA,it) + [6{t),Ys]dt = ^ 
[a,(t),n]=0, [A,(t),n]=0, [a*{t),Y,]=0, [5(t),n]-0, 

i.e. ctj{t), f3ij{t), aj{t), 6{t) £ At- This enables us to define a process Yt by 

dYt = (^ajit)dA*{t) - £*{V*aj{t))dt'^ + [p,,{t)dK,j(t) - £' {V* P,,{t)Vj)dt^ + 

(a*{t)dAj{t) - E\a*{t)Vj)dt^ , % = 0, 

i.e. we have the following splitting of Yt- 



(7.3) 



Yt=Yo+Yt+ f {£^{V*a,{s)) + £^ {V* P,,{s)V,) + {a*{s)V,) + 5{s))ds, (7.4) 



which in view of the following theorem is the semi-martingale splitting of Yt. The process Yt is called 
the innovating martingale of the observed process Yt . 

Theorem 7.3: The process {Yf}t>o is a martingale with respect to the titration {T,t}t>o of and 
the measure Pp, i.e. for all t > s > we have: E,p{Yt) = Yg. 

Proof. We need to prove that for alH > s > : Wp{Yt — Yg) = 0- This means we have to prove for 
ah t > s > and £■ e S^: 

E;iYt - Y,){uj)Yp{duj) = Q ^ [ {Yt- Ys){Lu)Fp{du) = ^ 

E JE 

Ep(YtE~Y,E- (^£''{V*ajir))E + f3,,{r)V,)E + {a* ir)V,)E + 6{r)Eyr^ = ^ 

p\YtE) - p^YsE) = ^ p'-(^£^{v;a,ir))E + (V* P,, ir)V,) E + {a* {r)Vj) E + 5(r)i;) dr. 

For t = s this is okay, so it remains to be shown that for alH > s > and _E G S^: 

dp\YtE) ^ p'(^£\V*aj{t))E + £\V*P,j{t)Vj)E + £'{a*{t)Vj)E + 6{t)Eyt 
dp"{Ut*YtEUt) = p*{£\V*a,{t))E + £\V* {t)Vj) E + £\a*{t)Vj)E + 6{t)E]dt. 
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We define: Zi{t) := U^, Z2{t) :— YtE and Z^lt) :— Ut then we find, using the notation below Theorem 
O dp°{U^YtEUt) = p°i[l] + [2] + [3] + [12] + [13] + [23] + [123]). Remember p° = 0®'=, i.e. we 
are only interested in the dt terms, since the vacuum kills all other terms. The terms [1], [3] and [13] 
together make up the usual Lindblad term and since L(l) = we do not have to consider them. 
Furthermore, term [2] contributes U^5{t)EUtdt, term [12] contributes UfV*aj{t)EUtdt, term [23] 
contributes a* {t)Vj EUfdt and term [123] contributes UfV* Pij{t)VjUtdt, therefore we get 

dp°{U:YtEUt) = p"{u:a*{t)VjEUt + U;V* (3^J(t)VJUt + U;v;aj{t)EUt + U;5{t)EUt)dt = 
p'{a*{t)VjE + V*P^j{t)Vj + V*aj{t)E + 5{t)E)dt = 

* {V*aj {t))E + f * {V*(3^j {t)Vj)E + £* (a* {t)Vj)E + 6{t)E'^ dt, 
proving the theorem. □ 

Remark. In the probability literature an adapted process which can be written as the sum of a 
martingale and a finite variation process is called a semimartingale |30j. The Theorems 17.21 and 17.31 
show that Mt and Yt are semimartingales. 

We now represent the martingale Mt from Definition l7.1l as an integral over the innovating martingale 

(cf. m) by 

dMt = VtdYt (7.5) 

for some stochastically integrable process rjf, which together with equation (|7.4|l provides the link 
between dMf and dYf. We are left with the problem of determining rjt, which we will carry out in the 
next section for the examples of Section [21 and |31 Here we just give the recipe for finding rjt- 

Recipe. Define for all integrable adapted processes bt and Ct a process Bt in Coo by 

dBt = btdYt + ctdt. (7.6) 

These processes form a dense subalgebra of Coo- Now determine rjt from the fact that £* leaves p* 
invariant jH], i.e. for all Bt 

p'{£\BtX))=p\BtX). 

From this it follows that for all Bt 

dp°{u;Bt{£\X) - X)Ut) = 0. (7.7) 

We evaluate the differential d(U^ Bt{£*'{X) — X)Ut) using the quantum Ito rules. Since p^ = p® 0®'^ 
we can restrict to the dt terms, since the others die on the vacuum. We will use the notation below 
Theorem O with Zi{t) = [/;, Z2{t) = Bt, Z^it) = £\X) - X and Z4^{t) = Ut- The following lemma 
simplifies the calculation considerably. 

Lemma 7.4: The sum of all terms in which Z2 is not differentiated has zero expectation: p"([l] + 
[3] + [4] + [13] + [14] + [34] + [134]) ^ 0. 

Proof. The dt terms of [3] are U; Bt£* {L{X))Utdt and -U; Btr]t{£\V*aj)+£\V* P,jVj)+£\a*Vj))Utdt. 
Using the fact that f* leaves p* invariant we see that the term Bt£* [L{X)')Utdt cancels against the 
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dt terms of [1], [4] and [14], which make up the Lindblad generator L with a minus sign. The other 
term of [3] is cancelled in expectation against the dt terms of [13], [34] and [134], since 

p°([13]) = p\BtVtV*aj)dt = p'{£\Btr,tV*a,))dt = {Btr,t£\V* aj))dt 

p\[M]) = p\BtT^ta*V,)dt - p'{£\BtT^ta*V,))dt = p' {Bt7u£\a*V,))dt 

p°([134]) = p\Btr,tV* l3,jVj)dt = p' {£\Btr,tV* P,jVj))dt = p' {BtVt£\v: (3^jVj))dt. 

□ 

Using equation (|7.3|) . the fact that leaves invariant and the module property, we find that the 
term [2] has expectation zero as well 

p°([2]) = {htdYt{£\X) - X)) = -p' {bt£'iV;a,{t) + a*{t)V, + V* P.jV,){£\X) - X)) dt = 

- p' {bt£*{V*aj{t) + a*{t)Vj + V*l3,jVj)£\£'{X) - X)) dt = 0. 

Thus, only the terms containing no B^ nor Ct can contribute non-trivially. This leads to an equation 
allowing us to obtain an expression for rjt by solving 

p"([12] + [23] + [24] + [123] + [124] + [234] + [1234]) = 0. (7.8) 

Although this can be carried out in full generality, we will provide the solution only for our main 
examples, the photon counting and homodyne detection experiments for a resonance fluorescence 
setup, in the next section. 



8 Examples 

We now return to the example considered in Section [2| We were considering a 2-level atom in inter- 
action with the electromagnetic field. The interaction was given by a cocycle Ut satisfying equation 
(|6.3(l . The observed process is the number operator in the side channel, i.e. Yt = Ass{t). Therefore 
dYt = dAssit) - £\V;Vs)dt. Recall now the notation Zi{t) = [/;, Z2(<) = Bt,Z3{t) = £\X) - X 
and Z4{t) = Ut, their differentials are given by 

du: = u: {KdA^t) - v^dA^it) - ^{-tH + v:-v,)dt} 

dBt = 6tdA,,(t) + {ct - ht£\v:V,))dt 
d{£\X) - X) = T]tdKss{t) + (£'{L{X)) - r]t£\v:Vs))dt 

dUt = J2 {^-'^Kit) - KdA^it) -liiH + v:V^)dt}Ut. 

Following the recipe of the previous section we now only have to determine the dt terms of [12] , [23] , [24] , 
[124], [123], [124] and [1234]. All of these terms are zero in expectation with respect to except for 
[124] and [1234] 

p°([124]) = p'{u:hV:{£\X) - X)VsUt)dt 
p"([1234]) = p\u:btiitV:VsUt)dt. 
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For all bt the sum of these terms has to be in expectation, i.e. 

ybt: p'(^bt{v:{£\X)-X)Vs+mV:V,^^dt = ^ 

ybf. p'(^£*(^bt(^v:{£\X)-X)Vs+vtV:V,^^'^dt^O ^ 

ybt: p'(^bt[£'iX)£\v:Vs)-£'{V:XV,)+vt£\v:Vs))^dt = ^ 

Substituting the expressions for 774 and Yt into equation H7.2|l we obtain the Belavkin equation for 
photon counting in the side channel 

d£\X) ^ E'{L{X))dt + - ^HX)) [dkUt) ~ £'{V:Vs)dt) . (8.1) 

Now recall that £\X) = pi{X,), i.e. it is the function Qt ^ C : uj ^ pI{X^). For &\\ X B ^ M2, 
the M2 valued function X, is the constant function lo ^ X. Therefore for all X in the Belavkin 
equation H8.1|l is equivalent to 

dp\{X) = pi{L{X))dt + ( ^'[[^^^y - P\{X)) (dAUt) - pliV:Vs)dt), 

which is equivalent to the Belavkin equation of Section 12 equation (|2.5() . In simulating the above 
equation we can take for Yt = Assit) the unique jump process with independent jumps and rate 
pi{V*Vs), since Ass{t) — /g pl{V*Vs)dr has to be a martingale. 

Let us now turn to the homodyne detection scheme which we already discussed in Section |3| The 
observed process is now Yt — ^^(i) = ^lift) + ^s(/t) (see Section 01 for the definition of ft). This 
means the innovating martingale Yt satisfies dYt = e"^'dA*(t) + e-*'**dAs(t) - £^e"l"V* + e~"^'Vs)dt, 
where (j)t = 00 + '-^lot with loio the frequency of the local oscillator. Therefore we find different 
differentials for Bt and £*{X) — X than we had in the photon counting case 

dBt = bt{e'^*dA:{t) + e~''^'dA,{t)) + (q - btS'ie'^'V; + e-'*'V,))dt 
d{£\X)-X) = T^t{e'*'dAl{t)+e-'^'dAs{t)) + (f*(i(X)) - Tjt£\e"^*V: + e-~'^'Vs))dt 

Following the recipe of the previous section we now only have to determine the dt terms of [12], [23], 
[24], [124], [123], [124] and [1234]. All of these terms are zero in expectation with respect to p", except 
for [12], [23] and [24] 

^ dt 



dt. 



p°([12]) = p''{u:^t'^V:bt{£\X) - X)Ut 

p''{m)^ p\u*tbtmUt)dt 

p°([24]) = p°{uM£\X) - X)e-'t'^VsUt) 
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For all bt the sum of these terms has to be in expectation, i.e. 

V6t: p'(^bt(^e''f'^V:{£'{X)-X) + {£\X)-X)e-'^'Vs+m)^dt = 4=^ 

ybt : p* (bt{- S^e^f^V^X + e-^f^XYs) + £*{e"*"V; + e-"'"Vs)£*{X) + f]t)yt = 



Vt 



Substituting the expressions for rjt and Yt into equation (|7.2|l we obtain the Belavkin equation for the 
homodyne detection scheme 

d£\X) = £^{L{X))dt + {£\e"l''V*X + e-^f^XYs) - £\e"'"V; + e^"^'Vs)£\X))x 
X {e"^'dA*{t)+e-"t"dA,{t)~£\e"'"V: +e-'''>"Vs)dt). ^^'^^ 

Now recall that £\X) = pi{X,), i.e. it is the function Vlt ^ C : uj ^ pI{Xu)- For all X e 6 = Ah, 
the M2 valued function X, is the constant function uj t-^ X. Therefore for all X in B, the Belavkin 
equation 18.2|l is equivalent to 

dpiiX) = pi{L{X))dt + {pi{e'^^V:X + e'^'^^XV,) - p\{e"^'V: + e-''^'F«)p* (X)) x 
X {e"'"dAl{t) + e~"l"dAs{t) - p\(e"^'V; + e~''^'V^,)di), 

which is equivalent to the Belavkin equation of Section [SJ equation H3.4|l . Since A*(/f) + ^^(/t) — 
/o pK^^'^^'V* + e~'^'^''Vs)dr is a martingale with variance t on the space of the Wiener process, it must 
be the Wiener process itself. 
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